The scalar glueball from a tadpole-improved action 
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The scalar glueball mass and the string tension are computed in lattice SU(3) gauge theory with the aim of 
establishing the effectiveness of the improved action approach in removing finite-spacing artifacts. 



1. INTRODUCTION 

The removal of lattice-spacing artifacts from 
numerical estimates of physical quantities is a 
central problem in lattice field theory. Taking the 
continuum limit by decreasing the coarseness of 
the mesh and extrapolating is a computationally- 
expensive task. The use of improved actions 
promises to be a much more efficient means of 
reducing cutoff contamination Q. The purpose 
of this work is to examine the effectiveness of the 
improved action approach in lattice SU(3) gauge 
theory by studying the low-lying glueball spec- 
trum. 

In this work, we present results for the scalar 
glueball mass m g and the string tension a using 
the tadpole-improved Luscher-Weisz action with 
lattice spacings a « 0.25 and 0.40 fm. Comparing 
to results from the standard Wilson action, we 
find a significant reduction in the finite-spacing 
errors in the ratio m g / y/a (setting the scale us- 
ing the string tension). The calculation exposes 
two problems with using the Liischer-Weisz ac- 
tion on a coarse lattice: first, the coarseness of 
the lattice in the temporal direction severely lim- 
its the number of statistically-useful correlator 
measurements and hampers the demonstration of 
plateaux in effective masses; secondly, the un- 
usual properties of the transfer matrix can re- 
duce the effectiveness of the variational method in 
diminishing excited-state contamination. These 
properties are due to the presence in the improved 
action of terms which couple fields separated by 
more than one time slice. 

First, the improved action is described in 
Sec. 0. Simulation details are then presented in 
Sec. H, followed by results and conclusions. 



2. THE IMPROVED GAUGE ACTION 

The perturbatively-improved QCD gauge ac- 
tion with tadpole improvement advocated in 
Ref. 0] was used in this study. The original per- 
turbative improvement was carried out long ago 
in Refs. J^H; the application of tadpole improve- 
ment H to this action is more recent. The action 
is given by 

S[U] = /3 pl 2 |Re Tr (1 - U pl ) 
pi 

rt 

+ A> g E i Re Tr 0- - u PS ), (i) 
pg 

where U p \ denotes the usual plaquette, L7 rt indi- 
cates the product of link variables about a planar 
2x1 rectangular loop, and U pg is the parallelo- 
gram loop. These terms are depicted in Fig. |l|. 

The coupling /3 p i is the only input parameter. 
The other couplings are computed in tadpole- 
improved perturbation theory by removing 0(a 2 ) 
errors in spectral quantities order-by-order in the 
QCD coupling. To one-loop order, these cou- 
plings are given in SU(3) by 

&t = -JK (1 + 0.4805 «□), (2) 

P pg = 0.03325 a n , (3) 

where the mean-field parameter uo and the QCD 
coupling an are given in terms of the measured 
expectation value of the plaquette as follows: 

u = (iReTrf/pi) 1 / 4 , (4) 
a a = -1.30362 lnu . (5) 
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Figure 1. Wilson loops in the improved action. 

Using identities from Ref. 0, the action S[U] was 
found in Ref. to be positive semi-definite for 
/3 P i > 6.8. 

Since the input couplings depend on the mea- 
sured quantity uo, the appropriate value for uq 
must be determined by tuning for each value of 
/3 p i. There are many ways this could be done. We 
chose to fix the /3 p i value and vary the input value 
of uo until agreement between the input and mea- 
sured values is reached. An easy and efficient way 
to do this is to use the Ferrenberg-Swendsen tech- 
nique || to simultaneously (in a single Markov 
chain) measure uo for a range of input values. 
These runs are very inexpensive. A sample of our 
tuning results follows: Uq(6.8) = 0.467, Uq(6.9) = 
0.480, t4(7.0) = 0.494, and ^(7.4) = 0.555. 

Due to the presence of terms which span two 
time slices in the action S[U] given in Eq. |l|, 
the concept of a transfer matrix is more compli- 
cated than for the simple Wilson action. Liischer 
and Weisz Q have shown that a transfer matrix 
can be defined, but it is in general not Hermi- 
tian; complex eigenvalues can occur, leading to 
damped oscillatory behaviour in correlation func- 
tions, and negative weight contributions can ap- 
pear in the spectral decomposition of two-point 
functions. Nevertheless, the mass spectrum can 
still be obtained from the exponential decay of 
correlators for large temporal separations. Care 
must be exercised, however, in applying varia- 
tional methods to extract masses from correlators 



at short time intervals: although invalid, strictly 
speaking, the variational method may still be 
practical whenever the coherence length is large 
enough such that the contributions from the prob- 
lematic high-energy states are negligible. 

3. SIMULATION DETAILS 

A set of 18 scalar glueball creation operators 
{Oi(t)} were constructed from combinations of 
Wilson loops on time-slice t with appropriate 
quantum numbers. Due to the coarseness of our 
lattice, fuzzing was not performed. We then ob- 
tained Monte Carlo estimates of the following ma- 
trix of correlators: 

C y (r) = 52{{Oi(t) Ojit+r)) - (Oi) (Oj)} . (6) 
t 

The simulations were run mainly on DEC-a and 
Hewlett-Packard workstations; 3000 CPU hours 
of Cray T3D time were also used. The T3D was 
employed essentially as a task farm consisting of 
64 isolated workstations. The gauge fields were 
updated using both Cabbibo-Marinari [CM] pseu- 
doheatbath and overrelaxation [OR] techniques. 
Updates with the improved action were found to 
be five times more costly than for the Wilson ac- 
tion. One CM and one OR sweep were performed 
between correlator measurements, and the results 
were combined into bins of 1000 to eliminate au- 
tocorrelations. Signals were enhanced by exploit- 
ing Euclidean lattice symmetries. Results were 
obtained using f3 p \ = 6.8 on a 6 4 lattice (334 bins) 
and p \ = 7.4 on an 8 4 lattice (174 bins) with pe- 
riodic boundary conditions. One simulation using 
the Wilson action with ftw = 5.5 was performed 
on an 8 4 lattice (104 bins). 

To determine m g , the variational method H 
was first used to select optimal combinations of 
glueball operators O op t = J2j c j Oj '■> this involved 
choosing the Cj 's such that C op t (t op t)/ C op t (0) was 
a maximum, where C op t is the correlator asso- 
ciated with O op t- The optimizations were per- 
formed using i opt — 1,2. The effective masses as- 
sociated with these optimal correlators were then 
examined. The coarseness of the lattice severely 
limited the number of time separations for which 
signals could be obtained; for this reason, a so- 
phisticated correlated analysis was not possible. 
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Figure 2. Effective mass for /3 p i = 6.8 



Figure 3. Effective mass for (5 P \ = 7.4 



Effective masses for glueballs in other channels 
were also studied, but no acceptable signals were 
obtained. 

In order to set the scale, we also measured the 
string tension a. This was done by extracting 
the static quark potential V(r) for various r from 
ratios of appropriate Wilson loops. The results 
were then fit to the form 

V(r)=ar-^- + b, (7) 

in order to obtain a 2 a. Only on-axis r vectors 
were used. 

4. RESULTS 

Our results for the on-axis string tension are 
a 2 cr(6.8) = 0.91(2) and a 2 er(7.4) = 0.36(1). Sys- 
tematic uncertainties associated with the form of 
V(r) used to extract a dominate over statistical 
uncertainties; the errors quoted are estimates of 
the systematic effects. Taking s/a = 440 MeV, 
the lattice spacing can be set, yielding a(6.8) = 
0.427(5) fm and o(7.4) = 0.269(4) fm. 

The effective glueball mass plots are shown in 
Figs. H and |L Circles (crosses) indicate results 
extracted from the correlator of the glueball op- 



erator obtained by optimizing on time-slice 1 (2). 
Note that the effective masses do not decrease 
monotonically; this behaviour is a result of the 
presence of the two-time-step interaction in the 
improved action. Evidence of a plateau is clearly 
visible in the (3 P \ = 7.4 plot; large errors shed 
some doubt on the existence of a plateau in the 
/? p i = 6.8 plot. It is interesting to note that the 
plateau used in Ref. [|| to obtain m g at f3w = 6.4 
lies in the range 0.2—0.6 fm, whereas t = 1—3 cor- 
responds to 0.3 — 0.8 fm for our (3 p \ = 7.4 results, 
and t = 1 — 2 corresponds to 0.4 — 0.9 fm for the 
/3 p i = 6.8 data. 

Finite- volume effects on our glueball masses are 
expected to be small. Using Liischer's formula || , 
we infer that our mass estimates lie within i% of 
their infinite volume limits. Note that our vol- 
umes 7(6.8) re (2.6 fm) 4 and 7(7.4) « (2.2 fm) 4 
are of similar size to those used in Ref. || . 

We are confident that the t = 2 effective mass 
value shown in Fig. |^ gives a reliable determina- 
tion of the glueball mass at f3 p \ = 7.4. This value 
is shown in Fig. |] where it is compared with glue- 
ball determinations from various earlier works us- 
ing the Wilson action. Our glueball mass using 
the Wilson action with P\y = 5.5 is also shown in 



5. CONCLUSION 
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Figure 4. Glueball mass against string ten- 
sion: Wilson and improved actions. The circle 
is our /3 p i = 7.4 result; the square indicates our 
(3w = 5.5 Wilson action result. The diamonds are 
Wilson action results from Refs. ju|^L2|,|J. Tnc 
straight line is a fit from Ref. |ll]] to the Wilson 
data. 



this plot and agrees with a previous determina- 
tion |l^]. Comparing the Wilson f3w = 5.5 result 
with the P p i = 7.4 ratio from the improved action, 
one observes a remarkable three-fold reduction in 
the finite-et errors. 

We are less confident that the t = 1 effective 
mass shown in Fig. reliably represents the glue- 
ball mass at f3 p \ — 6.8. However, if we use this 
value and convert to physical units using y/a = 
440 MeV, we obtain m s (6.8) = 1.38(7) GeV, com- 
pared with m s (7.4) = 1.31(11) GeV. Results at 
more (3 p \ values are needed before one can safely 
comment about scaling properties. Note that the 
above values are low in comparison to the a — ► 
Wilson results. 



We computed the scalar glueball mass and 
the string tension using the tadpole-improved 
Liischer-Weisz SU(3) action on two coarse lattices 
and found a significant reduction in the finite lat- 
tice spacing errors in the ratio m g /y/a. The large 
lattice spacings in the temporal direction limited 
the number of statistically-useful effective mass 
measurements. In the future, we intend to cir- 
cumvent this problem by working on asymmetric 
coarse lattices in which the spacing in the time di- 
rection is reduced. We expect that the powerlaw 
increase in simulation cost will be more than off- 
set by the exponential increase in signal-to-noise. 
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